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Abstract 



A nonperturbative renormalization group equation describes how the momentum depen- 
dent four-quark vertex depends on an infrared cutoff. We find a quasilocal one-particle 
irreducible piece generated by (anomaly-free) multi-gluon exchange. It becomes impor- 
tant at a cutoff scale where scalar and pseudoscalar meson-bound states are expected to 
play a role. This interaction remains subleading as compared to the effective one-gluon 
exchange contribution. The local instanton induced four-quark interaction becomes dom- 
inant at a scale around 800 MeV. In absence of a gluon mass the strong dependence of 
the one-gluon exchange on the transferred momentum indicates that the pointlike interac- 
tions of the Nambu-Jona-Lasinio model cannot give a very accurate description of QCD. 
A pointlike effective four-quark interaction becomes more realistic in case of spontaneous 
color symmetry breaking. 



1. Introduction 



The Nambu-Jona-Lasinio (NJL) model |3J has often been used as a simple model for low- 
momentum strong interactions. After inclusion of the chiral anomaly it gives a relatively 
successful description of the chiral symmetry aspects of QCD and the associated properties 
of mesons f2|. A renormalization-group treatment based on exact flow equations has 
provided a consistent picture of the chiral aspects of the high temperature phase transition 
for two light quark flavors || . On the other hand it is well known that crucial properties 
of long-distance QCD such as confinement cannot be accounted for by this model. The 
NJL-model shares the chiral symmetries of QCD whereas color symmetry appears only 
as a global symmetry. It is based on a pointlike four-quark interaction 

rff = / d 4 * {0 a Yd^t + (ft^) } , (i.i) 

where A NJL indicates that only fluctuations with momenta q 2 < A^ JL should be included 
in the functional integral. (Here i,j = 1,2,3 and a, b = 1,2,3 are the color and flavor 
indices, respectively.) The quartic interaction may be supplemented by a local anomaly 
term which is sixth order in the quark fields for three light flavors. Also extensions with 
four-quark interactions in the vector channel have been investigated p|. 

From the perspective of QCD, integrating out the gluon fields as well as all fermion 
fluctuations with (covariant) momenta larger than Anjl should lead to an effective action 
containing pieces like (|1.1|) , together with possibly quite complicated additional multi- 
fermion interactions. The question arises under which circumstances and for which prob- 



lems the NJL model ( |1 . 1| ) can be considered as a good approximation to the multi-quark 



interactions obtained from QCD. As a first step towards an answer we compute in this 



paper the four-quark vertex with the index structure appearing in Eq. (LI]). An inves- 
tigation of its momentum dependence serves as a test if the pointlike interaction of the 
NJL model is reasonable. 

Our method to explore the physics at the nonperturbative scale Anjl ~ (500 — 800) 
MeV is based on truncations of the exact nonperturbative renormalization group equa- 
tion for the effective average action H, H, H, 0, |§. An earlier study || has followed the 
running of the momentum-dependent effective four-quark interaction with initial con- 
ditions "mimicking" QCD. In contrast, we start in this paper at short distances from 
standard perturbative QCD. We observe that a one-particle-irreducible (1PI) four-quark 
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interaction is generated by box-type diagrams and becomes indeed approximately point- 
like for scales of the order A^jl. Nevertheless, the contribution of the one-gluon exchange 
to the effective four-quark vertex retains the characteristic momentum dependence of a 
particle exchange in the t-channel. This contribution dominates quantitatively over the 
NJL-interaction. Another pointlike interaction arises from instanton effects. It becomes 
dominant at a cutoff scale around 800 MeV. We show the different contributions to the 
four-quark vertex in Fig. 2, where we present the case of massless gluons. For massless 
gluons a pointlike NJL-interaction of the type ( |1 . 1|) seems not to be a realistic approxi- 
mation to the effective four-quark vertex in QCD. The hypothesis that gluons acquire a 
mass My m 850 MeV from spontaneous color symmetry breaking irD| is discussed later 
(Sect. 5). For massive gluons the approximation of pointlike four-quark interactions can 
be defended. Nevertheless, any quantitatively reliable approximation should retain the 
one-gluon exchange in addition to the instanton interaction and the NJL-interaction ( |1.1| ). 
This is the setting used for gluon- meson duality pTj| . 

The paper is organized as follows. In Sect. 2 we derive a truncated nonperturbative 
flow equation for the scalar-exchange-like effective momentum-dependent four-quark in- 
teraction in QCD. This effective one-particle irreducible interaction arises from box- type 
diagrams. The "box interaction" is relevant for the physics of scalar and pseudoscalar 
mesons and therefore for spontaneous chiral symmetry breaking. In the limit where it 
becomes independent of the momentum exchanged in the t-channel it can be described 
by meson exchange. In Sect. 3 we study, by numerical methods, the above-mentioned 
flow equation. We also present the perturbative expansion of the flow equation and its 
solution and comment about some checks on the validity of our approximations. In Sect. 
4 we compare the box interaction and the local instanton-induced four-quark interaction. 
In Sect. 5 we discuss the hypothesis of spontaneous color symmetry breaking, and a 
summary and conclusions are drawn in Sect. 6. 



2. Nonperturbative flow equations 

i) EFFECTIVE AVERAGE ACTION 

The effective average action Tk @] generalizes the generating functional for the IPI 
correlation functions F (effective action) for a situation where only fluctuations with (co- 
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variant) momenta q 2 > k 2 are included in the functional integral. It obtains by introducing 
an infrared cutoff ~ k into the original functional integral such that for k — > one re- 
covers the effective action T = Tj,^ . For a given at nonzero k the computation of T 
from Tfc involves only the fluctuations with q 2 < k 2 . This is the setting one needs for the 
effective NJL model. Our aim is therefore the computation of Tj. for k = A^jl- Since 
we do not know the precise value of Anjl we will consider arbitrary k in the appropriate 
nonperturbative range. More precisely, the NJL model is formulated in terms of quark 
fields. Since we start with the QCD action involving both quarks and gluons, the gluons 
should be integrated out. This may be done by retaining the infrared cutoff only for the 
fermionic fields. One then has to solve the gluon field equations obtained from as a 
functional of the quark fields and reinsert the solution into T^. [IT, T2\. Typically, we will 



not integrate out the gluon fluctuations with momenta smaller than k in this paper. Our 
results will nevertheless be a good guide for the qualitative characteristics of the effective 
four-quark interactions. 

The effective average action Tk[A, c, c, if), if)] for gluons A, ghosts c and quarks if) obtains 
by a Legendre transform from a functional integral involving the QCD classical action and 
an infrared regulator: 

S k [A,C,C,lf),lf)] = S gauge [A, C, c] + AS k ,gauge[A,C,c] + S F [A,ip,ip] + AS kF [ip,ij] . (2.1) 

Here the fermionic action S F is given by 

S F [A^3] = J d 4 x rfjfip^A])^ , (2.2) 

with (D^A})^ = Sjd^ —igA z (T z ') i 3 the covariant derivative in the fundamental represen- 
tation. Here a denotes the flavor and we use two sorts of color indices z — 1, . . . , N 2 — 1 
and i,j = 1, . . . , N c for a gauge group SU(N C ). We consider, for simplicity, Nf massless 
quarks: a = 1, . . . , Nf. Throughout this paper we work in four- dimensional Euclidean 
space, so that the matrices 7 M in the previous equation are the (hermitean) Euclidean 
7-matrices. The gauge part S gauge consists of three pieces 

S g auge[A, c, c] = S YM [A] + S gf [A] + S gh [A } c, c] , (2.3) 

where Sym is the Yang-Mills action, given by: 

S YM [A] = yd i xF* u F* u , (2.4) 
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with 

F; v = d,K - d v A; + gf"° A vAZ . 
The gauge-fixing part S g f reads 

S9f[A] = ~k J ' d * x d » A 'A A i > 

where a is the gauge parameter, and the ghost action S g h is given by 



(2.5) 



(2.6) 



S gh [A,c,c] = J d 4 xd,c z (5 zw d, + gr w A 



(2.7) 



The infrared cutoff modifies the propagators. For the fermions, the explicit expression 
ASkF[ip, V>] reads in momentum space || O 



A5 JfcF [V,V] 



d 4 g - 



(2tt)< 



(2.8) 



where we choose here 



e -g 2 A 2 



1 _ e -q 2 /k 2 

In particular, the (regularized) free fermion propagator turns out to be: 



(2.9) 



G?2{q) 



ab 



5 ab 5, 



7% 



g 2 [l + r fcF (g2)] 



5 ab 5^ 



l-e^ 2 / fc2 



(2.10) 



The gauge infrared cutoff term ASfc ifl£mfle [,<4., c, c] is the sum of a gluon infrared cutoff term 
and a ghost infrared cutoff term: ASk, g auge[A, c, c] = ASfc^lyl] + AS^^fc, c]. The explicit 
expression of the gluon infrared cutoff term AS^^A] is given by || |[ 



where 



AS kA [A\ 



d 4 p 
2 J (27r? 



A*(-p)i#(PK(p) 



<v + 



- 1 



a: 



^2 



^a(p 2 ) 



(2.11) 



(2.12) 
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with 



Rk A (p 2 ) = p 2 



-p 2 /k 2 



1 _ e -p 2 /k 2 

Finally, the ghost infrared cutoff term ASk, g h[c,c] reads 



AW,*] = / 7S1 c z (-p)R k , gh (p 2 )c z (p) 



(2vr 



where 



- P 2 /k 2 



R k , gh (p)=P T _ e _ p2/k2 . 
With this choice, the regularized ghost propagator becomes 



Z.V 



V 2 + ^(p 2 ) 



1 _ e -P 2 A- 2 
P 2 



(2.13) 



(2.14) 



;2.15) 



(2.16) 



It is well known || 0, || that the presence of the infrared cutoff explicitly breaks gauge 
and hence BRS invariance: in particular, a gluon mass term appears. A background field 
identity |J or modified Slavnov- Taylor identities || [13| can be derived, which are valid for 
k > 0. They guarantee the BRS invariance of the theory for k — > 0. Within a perturbative 
expansion of T k) one can derive an equation for a fc-dependent mass term for the gauge 
fields implied by the modified Slavnov- Taylor identities. For pure Yang-Mills theories 



(i.e., with no fermions), regularized with the above infrared cutoff term (|2.11|) - (|2.13|) , one 
finds the following result, at the leading perturbative order 
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rn 



kA 



1287T 2 



k 2 {a - 1) . 



(2.17) 



Here a is the gauge parameter, introduced in Eq. (|2.6j), and g k is the renormalized 
(/c-dependent) coupling constant. In the theory with regularized fermions that we are 
considering, a contribution to the gluon mass also appears from the fermion part of the 
gluon self-energy. One finds, at the leading perturbative order: 



m 



kF 



W f g 2 k 

327T 2 



(2.1* 
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When including the gluon mass, the regularized gluon propagator becomes, in the Feyn- 
man gauge a — 1: 



p 2 + m 2 k + R kA (p 2 



where 



m l = m kA + m kF 



(2.19) 



(2.20) 



(Indeed, m\ A = at the perturbative order 0(g 2 ) in the Feynman gauge.) In other words, 
using the expression ( p,13| ) for R k A{p 2 )'- 



-HI 



1 _ e -P 2 A 2 



p 2 + m|(l — e _p2 / fc2 ) 



X21) 



We shall concentrate on the part T k F \ijj,ijj] of the effective average action which de- 
pends only on the fermion fields: 



T k [A, c, c, ^, ^ = rf } [V, $ + T k A > [A, c,c,j>,tf\. 
In general, it can be expanded |TT], [TIJ as 

rf ) = rg + rS + . . . , 



(A)r 



;2.22) 



(2.23) 



where r^^fV'jV'] contains terms with 2n fermion fields, i.e., of the form ~ (ipif)) n . For 
n — 1, we have the bilinear termF] 



d 4 g 



(2.24) 



Similar to Refs. |9|, |12| , we truncate the effective action by omitting T^/ k [ip, ip] with n > 3 

' ( F} r n 

and we concentrate on the behavior of the four-quark interaction (n = 2). For T 4k [ip,ilj] 
we adopt a chirally-invariant parametrization |[T2|| : 

X {A ct (Pi,P 2 ,P3,P4).Ht(P1,P2,P3,P4) + A,(pi, p 2 , p 3 , Pl)M p {p 1 , P2, P3, Pi) 

+*p(Pi,P2,P3,Pa)M p (pi,P2,P3,P4) + ...} , (2.25) 



*We omit in this paper the running of the fermion wave function renormalization. 
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where the four-fermion couplings A; have the property: 



Ai(pi,p 2 ,P3,P4) = Ai(-P4,-P3,-P2,-Pl) , 

and the usual invariant quantities s and t read 

S = (pi + p 2 ) 2 = (P3 + P4) 2 , t = (Pi - P3) 2 = (P2 - P4) 2 

The four-fermion operators M.i are defined as: 



£.26) 



£.27) 





1 

~2 


Wa 


(~Pl)^(p 2 )] ^ 








1 

+ 2 


?a 


(-Pi)7 5 ^-(p 2 ) 




^(p 4 )7 5 ^(-P3)^ 




M p (p 1 ,p 2 ,P3,P4) 


1 

= + 4 


?a 


(-PiW-(p 2 ) 




^(P4)7^ a (-P 3 )' 






1 

+ 4 


?a 


(-pi)7VV>}(Pa) 




^(P4)7m7 5 ^(- 


-Ps) 


M p (p U p 2 iP3iP4) 


1 r 

2iV c i 


C(-Pi)7^ a ( 


-Ps) 


[Vfeh/^fe) 



(2.28) 

If A CT does not depend on t, the term .M CT describes the exchange of color-singlet, flavor- 
nonsinglet spin-0 scalar quark-antiquark states in the s channel: they correspond to the 
spin-0 pseudoscalar and scalar mesons. Similarly, A4 P describes the exchange of color- 
singlet, flavor- nonsinglet spin-1 vector states in the s channel: they correspond to the 
spin-1 vector and axial vector mesons. Finally, M. p describes the propagation of a color- 
and flavor- singlet spin-1 vector state in the t channel: these are the quantum numbers 
of the pomeron. We concentrate on the analysis (both analytical and numerical) of the 
evolution equation for the scalar-like effective four-fermion coupling A CT , which multiplies 
A4 a and is relevant for the physics of scalar and pseudoscalar mesons and therefore for 
spontaneous chiral symmetry breaking. The dots in Eq. ( |2.25| ) stand for more complicated 
momentum structures, the contribution of the chiral anomaly (to be discussed in Sect. 4) 
and terms breaking the chiral flavor symmetry. 

ii) FLOW EQUATIONS 

The exact flow equation for can be written as a formal derivative of a renormaliza- 
tion group improved one-loop expression H 



dk^k 



1. 



-STr 



4in(rl 2) + J R fc ; 



(2.29) 
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(2) ~ 

where T k is the second functional derivative with respect to the fields and dk acts only 

(2) 

on Rk and not on T k , i.e.: 

4 = "aFa^ • (2 ' 30) 

Here Rk is to be interpreted as a matrix Rk = diag (-R^a, Rk, g h, Ruf)- The trace in- 
volves a summation over indices and a momentum integration and "STr" indicates an 
additional minus sign for the fermionic part. Flow equations for n-point functions follow 
from appropriate functional derivatives of Eq. ( |2.29| ) with respect to the fields. For their 
derivation it is sufficient to evaluate the corresponding one-loop expressions in presence 
of the infrared cutoff (with the vertices and propagators derived from Tk) and then to 
take a formal dk derivative. This permits the use of (one-loop) Feynman diagrams and 
standard perturbative techniques. 

Using these techniques, we have derived the flow equation for the four-fermion (1PI) 
coupling Ao- by computing the relevant one-loop diagrams from the effective average action 
Tk[A,c,c,ip,^IJ\: these diagrams are shown in Fig. 1. Our truncation includes only two 
vertices: i) the usual quark-gluon vertex and ii) the four-fermion (1PI) vertex \ a M. a . We 
have extracted from the relevant diagrams the contributions to the cr-type four-fermion 
vertex M. a and then performed a formal dk derivative in the sense of Eq. ( |2.30| ). In this 



picture (which is different from the one used in Ref. ||) the gluons are not integrated 
out. Therefore, in addition to the contributions to the flow equation coming from the 
infrared cutoff in the fermion lines inside a given Feynman diagram, one has to consider 
also the ones from the infrared cutoff in the gluon lines. 

iii) TRUNCATIONS AND THE FLOW EQUATION FOR A CT 

In general, the function \ a can depend on six independent Lorentz-invariant products 
of the momenta A CT = A CT (s, t, p 2 , p\-, P35 pi) • We have derived the evolution equation of 
Ao- for the particular configuration of momenta (^1,^2,^3,^4) = (p, — P, — P,p), so that 
s — 0, t — Ap 2 . In order to simplify the momentum structure, we have made a further 
truncation, namely that A CT depends only on the invariant t: X a ~ A CT (t). Furthermore, 
our truncation restricts the quark-gluon interactions to minimal coupling with a coupling 
constant depending on the gluon momentum. We use the classical form of the (regularized) 
gluon propagator[[ Our truncation omits the influence of all multiquark interactions 
except those proportional A4 a ( |2.28|) . Moreover, we neglect the possible effects from a q 2 - 



^See [[l3| for an investigation of the g 2 -dependence of the gluon propagator by means of nonper- 
turbative flow equations. 
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and /c-dependent fermion wave-function renormalization Z^. 
With these approximations, we find the following result: 



kd k X a (t) 



l6N r . 



dl « ^ 1 



k 2 J (2vr) 4 [ V 2N c, 
xq 2 F k (q)G k (p - q) G k {j> + q) + F k (q)A k (p + q) 

1 



-ira k ({p + q) 2 )q 2 



2Nr 



A CT (V) e-o/ k F k {2p + q)G k {p + q) 



+e -(2 P+q) y k i Fkiq)Gkip + g) + Fkiq)Fki2 p + q)A k (p + q) 
(~ 3 



2N 2 



K((P - q?)F k {q) 2e-« ' k G k (p + q) + F k (q)A k (p + q) 



+K{{p - qf)K{{p + g) 2 )gV^/ fc2 F fc (g) . 



(2.31) 



Here we use the following shorthands: 
F k {q) 1 
G k (q) 

Mq) 



and employ 



1 - e-« 2 / fc2 



q 2 [l + r kF (q 2 )] 
1 



l_e~« 2 /fc 2 



q 2 + m 2 + R kA (q 2 ) q 2 + m 2 (1 - e"« 2 / fc2 ) 
[g 2 + ml(l-e-'? 2 / fc2 )] 2 ' 



d k F k (q) 
d k G k (q) 



dr kF (q 2 ) dF k (q) 
dk dr kF (q 2 ) 

dR kA (g 2 ) dG k (q) 
dk dR kA (q 2 ) 



__ P -9 2 A 2 



(2.32) 



(2.33) 



The first two terms in the r.h.s. of Eq. ( |2.31| ) (i.e., those proportional to ~ a\ or ~ a k \ a ) 
reflect the exchange of gauge bosons [see Fig. la) -j- lc)]. The last term (i.e., the one 
proportional to ~ A 2 ) corresponds to the result of Ref. || for the truncation A CT = \ a {t) 
[see Fig. Id)]. 
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iv) RG-IMPROVEMENT OF THE GAUGE COUPLING 

The flow equation ( |2.31| ) is "RG-improved" in the sense that we use everywhere a 
"running" coupling constant oik(q 2 ) = g 2 (q 2 )/An. Instead of computing ak(q 2 ) by its own 
flow equation, we use here an "ad hoc" relation to the two-loop running gauge coupling: 



0.(11) 



An0 o \n(^/A 2 QCD ) 



1 - 



ft ln[ln(/x 2 /A| CT )] 
(51 H^I^qcd) 



;2.34) 



Here 0q and 0\ are the two first coefficients of the QCD /5-function 



0(9) = » 



dg_ 



- Pig 5 + 



(2.35) 



which, for a gauge group SU (N c ) and Nf quark flavors, are given by: 



0o 



(An) 



11N C - 2N f 



0i 



(An) 



34 2 /l3A 2 -3\ Ar 

— N 2 - N 

3 c \ 3N r l ■ 



(2.36) 



such that, for N c 



3 and Nf = 3: 



a(n) 



4vr [ _ QA \n{\n(^ 2 /K 2 QCD )] 

VH^/Xqcd) [ 81 \n(^/A 2 QCD ) 



(2.37) 



For the relation between q 2 , k 2 and /x 2 , we use two different ansatze: (A) fi 2 = q 2 + k 2 , and 
(B) fi 2 = k 2 . The first one (A) is motivated by the fact that both q 2 and k 2 act as effective 
infrared cutoffs such that essentially the larger one counts. The second one (B) neglects 
the dependence of the vertex on the gluon momentum. The overall mass scale needs in 
addition the unknown ratio between Aqcd i n the scheme relevant for the flow equations 
(ERGE-scheme) to the MS scheme. We therefore report all quantities with dimension 



of [mass] n in units of (GeV • Aqcd/ '^Sf)"; "where A^ is the fundamental mass scale of 



(3) 



MS 



MS 



QCD in the MS renormalization scheme (two-loop), with Nf = 3. The experimental value 
a s (Mz) = 0.118(3) corresponds to A^- ~ 360 MeV. For comparison we will also use a 



fixed value of a corresponding to the scale \i = k c = 1.5 • Aqc^/A^- GeV, i.e., ~ 0.34 
[ansatz (C)]. 
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ONE-GLUON-EXCHANGE CONTRIBUTION 



The term X a Ai a describes only the contribution to the effective cx-type vertex coming 
from the 1PI graphs. A four-quark interaction with a similar structure arises also from 
one-gluon exchange which is not one-particle irreducible and therefore does not contribute^ 
to A CT in the present formulation. Indeed, the four-fermion vertex induced by the one-gluon 
exchange is given by [|T2" 



Ar (l9) 



x2na k ((p 1 -p 3 ) 2 ) ■ G k (p x - p 3 ) ■ M(p u p 2 ,P3,P4) , 



(2.38) 



where 



Ij.bl 



(2.39) 



By an appropriate Fierz transformation, one can show that 



M = M a + M p + M 



p • 



(2.40) 



where -M^, M. p and M. v have been defined in Eq. ( 2.28 ). The total contribution in the 
cr-channel is therefore ~ (A CT + with 



X^\t) = 27ra k (t) 



1 - e-*/ fc2 



t + m 2 k (l - e-^ 2 ) ' 



(2.41) 



This contribution remains nonlocal for k — > if m k tends to zero in this limit. The size 



of the 1PI contribution A CT should be compared with the effective coupling Q2.41|) . 



*In the formulation || |TT|, g the one-gluon exchange part is included in A 
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3. Scale dependence of the four-quark interaction 



i) PERTURBATION THEORY 

In perturbation theory the 1PI four-quark interactions arise from "box" and "cross" 
two-gluon exchange diagrams [see Fig. la) and lb)]. For the part with the structure M. a 
we find the perturbative value (t = 4p 2 ) 

AW(f) = lQn 2 N c (l - 2^2) / ^j-MiP - + ^ 

xq 2 [F k (q)] 2 G k (p-q)G k (p + q) . (3.1) 

In particular, for t = and momentum-independent a k [ansatz (B)], the momentum 
integration can be performed explicitly. One finds, for N c = 3 and in the limit of a 
vanishing gluon mass term m 2 . = 0: 

\<*\t = 0) = ^-K a . (3.2) 

Here the constant 

d 4 q 



K„ = StT 8 * 2 / ' ^q 2 [F k {q)f[Gl{q)f , 



(27T) 

o 1 l-e-" 2 / fe2 

k q 2 + R kA {q 2 ) <T 



depends on the precise formQ of the infrared cutoff. For our choice [see Eq. fl2.13|) ] one 
finds 



K a — 2 • [5 In 2 — 3 In 3] — 0.34 . (3.4) 
For the purpose of comparison with our results in the following figures, we quote the value 

MS 



corresponding to k = k c = 1.5 • Aqc^/A^. GeV (a kc ~ 0.34): 



A^(t = 0)| fc=fcc ~ 8.7 x 10- 2 (GeV ■ Aqcd/A^Y ■ (3.5) 



In order to get a feeling for the relative strength of the a-channel as compared to other 
types of 1PI four-quark interactions, we also present here the full one- loop 1PI interaction 
*For the fermion cutoff used in H one has to replace F 2 — > F/q 2 . 
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for zero external momenta (i.e., s — 0, t = 0). For an arbitrary gauge parameter a one 
obtains the effective (local) four-quark interaction 



AF)(p) 9k K 

4 < fc 64tt 2 A; 2 a 

9 9 3 2 
— I — a H — a 
4 2 4 



-\ (feVVtf) (^7m7 5 ^-) ~ (%tV V>?) (^7m7 5 ^) } • (3.6) 

We note that the first term has precisely the structure of the one-gluon-exchange diagram 
~ M. (|2.39| ). We also emphasize that there are several ways how to Fierz-transform the 



expression (|3.6| ). Our truncation of keeping only A,,- corresponds more precisely to state 
which invariants are neglected. 

It is straightforward to verify that the first term in the flow equation ( [2.31|) for A^ 



precisely corresponds to the formal d k derivative of the perturbative value (|3.1|). For t = 
the contribution to kd k \ a obtains by replacing in the r.h.s. of Eq. (|3.2|) the constant K a 
by -1K a ~ -0.68 . 

ii) INITIAL CONDITION FOR THE FLOW EQUATION 

In order to solve explicitly the flow equation ( j2.31j) , we need to fix the initial condition, 
i.e., the value of \ a (t) at a certain "initial" value k = k c . We assume that at the scale 
k — k c — 1.5 • Aqcd/A^ GeV, where perturbation theory is assumed to work well, the 
value of \ a {t)\k=k c is well approximated by the leading perturbative contribution. In other 
words, we use as initial condition for the flow equation ( p.31| ): 

X a (t)\ k=kc = X^(t)\ k=kc . (3.7) 



We have then solved numerically the flow equation (|2.31| ) with the initial condition ( j37T 



iii) RESULTS 

In Fig. 2 we plot the resulting A CT (t) (for N c = 3 and Nf = 3) as a function of k, for 
the particular case t = 0. The continuous line corresponds to a running gauge coupling 
according to the ansatz (A). The dotted line in the same figure represents X a (t = 0) as a 
function of k with the initial condition \ a (t)\ k=kc = 0, instead of (|3.1D, (|3.7|). We observe 
the insensitivity with respect to the precise choice of the initial condition. The dot-dashed 
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line is the value of X a (t = 0) as a function of k with the initial condition (|3.1|) , ( 3.7|) , but 
using everywhere (both in the flow equation and in the initial condition) the momentum- 
independent coupling constant defined by Eq. ( |2.37| ) with fi 2 = k 2 [ansatz (B)]. We 
observe a strong increase of X a at a scale k ~ 400 — 600 MeV. 

Comparison with perturbation theory is made in Fig. 3. We plot the value of X^(t = 
0), defined by Eq. ( |3.1| ), as a function of k (dotted line). The dashed line represents the 
value of \<f\t = 0), defined by Eq. Q, function of k, but using the momentum- 
independent coupling constant defined by Eq. ( |2.37|) with /i 2 = k 2 [ansatz (B)]. In this 
case, the integration in Eq. ( |3.1| ) can be performed explicitly (neglecting also the gluon 
mass rrik appearing in the function Gk) and one obtains (for N c = 3) the result fl3.2|) -( [3~4"l) . 
Finally, the dot- dashed line is obtained by using this expression for X<f\t = 0), but using 
the fixed coupling constant = ctk=k c — 0.34 [ansatz (C)]. 

In Fig. 4 we plot X a (t) (continuous line) as a function of t at the fixed value k = 0.42 • 
Aqcd/A^t GeV. At this scale the approximation of the box-type four-quark interaction 
by a local NJL-vertex ( |1.1| ) is reasonable. 

iv) CHECKING THE APPROXIMATIONS 

When inserting the effective action (|2.22|) , Q2.25|) into the flow equation (|2.29|) and 
extracting only the contributions to the cx-type four-fermion vertex Ai a , one finds that 
in the flow equation ( [2.31| ) for the effective coupling X a also other pieces proportional to 
~ XiXj appear: these additional pieces are generated by the four-fermion structures A4 P 
and Ai p and they are put to zero in our truncation. Since only the retained contribution 
~ X 2 is proportional to the color factor N c , while all other terms of the form ~ XiXj are 
suppressed, our truncated evolution equation becomes exact in the leading order in a 1/N C 
expansion. Moreover, we have checked that a multiplication of the term proportional to 
X 2 in the flow equation ( |2.31|) by a factor of ten does not strongly influence the results. 
The solution X a (t) of the modified flow equation is very near to the real solution of Eq. 
( |2.31| ) reported in Figs. 2, 3 and 4. In other words, the terms of the form ~ XiXj do 
not seem to play a fundamental role in the flow equation at the scales of k reported here. 
Neglecting the contributions from the p and pomeron channels does therefore not seem 
to be a too brutal approximation. 

We have also tried to test the approximation involved in using the two-loop expression 
for the running coupling constant. In order to do this, we have considered a different 
version of the running coupling constant in the flow equation ( |2.31 ): we have taken for 
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gl(p 2 ) the usual two-loop expression when gl(p 2 )\2-ioo P < 10, while gl(p 2 ) = 10 when 
9k(P 2 )\2~ioo P > 10. In other words, the new g\{p 2 ) stops running when it reaches the value 
g\ = 10. Also in this case, the solution X a (t) of the modified flow equation does not show 
any dramatic difference when compared to the solution of Eq. fl2.31| ) reported in Figs. 2, 
3 and 4. 

Finally, since the one-particle irreducible four-point function remains small in the 
momentum range of interest, a neglection of the higher-order anomaly-free effective mul- 
tiquark vertices seems well motivated. 



4. Comparison with the instanton contribution 



In Eq. ( |2.25| ) we have chosen to parametrize ^^2[ip, tf\ by using U(Nf) ® U(Nf) chirally 
symmetric operators Ai a , M- P and M. p . In this section we supplement the effect of the 
chiral anomaly. We add to the effective four-quark interactions an anomaly-induced term 



j~5| , which is ST/ (2) ® SU(2) ® U(l)v invariant, but not U(1)a invariant. 



In particular, for a nonvanishing mass m s of the strange quark, the instanton-induced 
fermion effective action for Nf = 3 includes a four-fermion contribution of the form [HJ 



r 



(4) 



4„ \(N c ) n (N c ) , x 
inst I yj. tb /\i ns t ^ 4 instK^ ) 



(4.1) 



where 



O 



(No) 
4,inst 



#4 - ±51% ) x 



x { [tf Uj ] [d k di] - [tfdj] \d k m] + [u^Uj] [d k -i b di] - [u'^dj] \d k j 5 Ul ) } , (4.2) 



This operator comes from ALg^ + AL^\ where I = instanton and A = anti-instanton, 
and the effective coupling constant involves an integral over instanton sizes p: 



A, 



inst 



6 



0.35 



k 3 Jo 



a(/i) 



2tt 



a(k/z) 



cxp 



2tt 



a(k/z) 



(4.3) 
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Here p—1 GeV is the renormalization scale for the fermion fields and [|T(| 17]: 



N c 



2tt 



a(l/p) 



2N C 



cxp 



2tt 



a(l/p) 



(4.4) 



with 



4.6exp(-1.68iV c , 
Cn < ~ **(N e -l)\(N e -2)\ ■ (45) 

For a(l/p) we take the (two-loop) running coupling constant ( |2.34j ) with p = 1/p. For 
the mass of the strange quark we have used the value m s ~ 150 MeV. For N c = 3 and 
Nf = 3 one has C3 ~ 1.51 x 10~ 3 and we use Q2.37| ). We have integrated the instanton 
length scale p between and 1/k, since k is an infrared cutoff momentum scale. The 
long- dashed line in Fig. 2 is the value of \ ins t (which is t-independent!) as a function of k. 
One concludes that the instanton-induced four-quark interaction becomes dominant at a 
scale k around (somewhat below) 800 MeV. 



5. Spontaneous color symmetry breaking? 



It has been suggested [10] that confinement can be described by the Higgs phenomenon. 
In this picture the expectation value of a color-octet quark-antiquark condensate Xo gives 
a mass to the gluons 

M 2 = g 2 effX 2 . (5.1) 

In the limit of three massless quarks the mass of all eight gluons is equal. Due to spon- 
taneous color symmetry breaking the physical electric charge of the gluons is integer and 
they carry precisely the quantum numbers of the known vector mesonsQp, K*,u. For this 
reason M is identified with the average mass of the vector meson octet My w 850 MeV. 
Confinement arises in this picture in complete analogy to magnetic flux tubes in super- 
conductors. A nonlinear reformulation of the Higgs phenomenon preserves manifestly the 



gauge symmetry and shows the equivalence of the confinement and Higgs pictures ]T8 
Instanton-induced six-quark interactions destabilize the "color symmetric vacuum" and 
have been proposed as the origin of dynamical color symmetry breaking. 

*See pf for a discussion of the singlet vector meson state. 
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The results in the preceeding sections correspond to a vanishing octet condensate 
Xo = 0. In order to get some information on the Xir dependence of the effective four-quark 
interactions we present in this section results for a nonvanishing value of the induced 
gluon mass term M 2 . Since the precise relation between g e ff in Eq. ( |5.1| ) and the running 
gauge coupling in our approximation (|2.37|) is nontrivial |10[ , we present results for two 
different settings. 



For the first setting (A) we simply add a constant gluon mass term M 2 = M v 



(850 MeV) 2 in the propagator such that: 

m 2 k = M v + m 2 kA + m 2 kF . (5.2) 

Moreover, we use a momentum-dependent running gauge coupling a(/i), with fi 2 = q 2 + 
k 2 + M v . For the second setting (B) we work at constant Xo = 140 MeV with a running 
perturbative gauge coupling Q2.34| ) g 2 = Ana, i.e., M 2 = g\x\- This replaces the gluon 
mass term in Eq. ( |2.20| ) by 

2 2 2,2,2 / rr o\ 

m k = g k Xo + m kA + m kF . (5.3) 



In this case, we use a momentum-dependent running gauge coupling a(/i), with /i 2 = 
q 2 + k 2 + M 2 k , where ~M\ = 47ra(/I = 600 MeV) • x\- 

In Figs. 5 and 6 we compare \ a and for different running gluon mass terms m| 
given by ( p.2| ) (A), ( |5.3| ) (B), and ( |2.20| ) (C). We observe that the gluon mass damps 
the increase of both A CT and A^ 19 -'. On the other hand, a pointlike approximation to 
the anomaly- free four-quark interaction becomes now better justified, especially for (A). 
Thanks to the additional infrared cutoff M 2 the flow can formally be followed towards 
k — » without encountering a problem from a diverging gauge coupling: the growth of 
the gauge coupling is reduced by the presence of the (nonvanishing) gluon mass term. 
One should not expect, however, that our results remain valid for very small k. Indeed, 
it has been argued that the quark wave function renormalization Z^ ik multiplying the 
quark kinetic terrr^] Q2.24p drops by a factor of about three as k reaches a few hundred 
MeV. This reflects the binding of three quarks into a baryon. As a result, the running 
gauge coupling g k increases by a factor of three and the renormalized four-quark couplings 
A CT and \£ 9 ^ become multiplied by a factor of nine. The running of the wave function 
^This wave function has been set — 1 for the present work. 
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renormalization becomes therefore important for low k. We have not included it here 
(except for its contribution to the running gauge coupling), and this is the reason why 
we do not explore very low values of k even though the infrared problems are formally 
absent in presence of a nonzero gluon mass. 

We finally note that, in presence of six quark interactions reflecting the axial anomaly, 
an octet condensate Xo 7^ will also result in additional effective four-quark interactions 
linked to the anomaly. They are not included in the present investigation. We esti- 
mate their order of magnitude to be comparable to (or even larger than) the instanton 
contribution discussed in Sect. 4. 



6. Summary and conclusions 

In this paper we have explored the effective Nambu-Jona-Lasinio-type four-quark inter- 
actions at the nonperturbative scale Anjl ~ (500 — 800) MeV. For this purpose we have 
used a method based on truncations of the exact nonperturbative renormalization group 
equation for the effective average action @, || § . More precisely, we have derived a trun- 
cated nonperturbative flow equation for the scalar-like effective momentum-dependent 
four-quark interaction in QCD. This type of interaction is relevant for the physics of 
scalar and pseudoscalar mesons and therefore for spontaneous chiral symmetry breaking. 

The nonperturbative renormalization group equation describes how the momentum 
dependent four-quark vertex depends on an infrared cutoff k. The initial value of the 
vertex at short distances is computed from perturbative QCD. We have studied the above- 
mentioned flow equation both analytically, in the perturbative regime, and numerically, 
in the nonperturbative regime. For the nonperturbative regime we have derived the 
dependence of the scalar-like four-quark interaction on the infrared cutoff scale k and on 
the momentum exhanged in the t-channel. The results are summarized in Figs. 2, 3 and 
4. 

We find that a quasilocal one-particle-irreducible (1PI) four-quark interaction is gen- 
erated and becomes indeed approximately pointlike for scales of the order Anjl, where 
scalar and pseudoscalar meson-bound states are expected to play a role. The "quasilocal- 
ity" of this "box-interaction" is put in evidence in Fig. 4, where the dependence on the 



19 



momentum- scale is shown. 

As discussed in Sect. 4, another pointlike interaction arises from instanton effects. In 
Fig. 2 a comparison is made between the strength of the box interaction and that of the 
local instanton-induced four-quark interaction. At scales of the order Anjl, the instanton 
interaction dominates. 

Nevertheless, the contribution of the one-gluon exchange to the effective four-quark 
vertex retains the characteristic momentum dependence of a particle exchange in the t- 
channel. We show the different contributions to the four-quark vertex in Figs. 2, 3 and 
4. The effective one-gluon exchange contribution remains substantially larger than the 



pointlike box- type interactions. In consequence, the pointlike interactions (1.1) of the 
Nambu-Jona-Lasinio model cannot give a very accurate description of QCD if the gluons 
are massless. Models based on a pointlike instanton induced four-quark interaction look 
more favorable. Still, the one-gluon exchange is not negligible. 

The results shown in Figs. 2, 3 and 4 have been obtained with the assumption of 
massless gluons (in the limit of vanishing infrared cutoff scale, k — > 0). In Sect. 5 we 
have discussed the intriguing hypothesis that gluons acquire a mass My « 850 MeV from 
spontaneous color symmetry breaking [jTU| . Our results for a nonvanishing gluon mass 



are shown in Figs. 5 and 6. In particular, from Fig. 6, where the dependence on the 
momentum scale is shown, it is clear that for this scenario pointlike four-quark interac- 
tions can be defended as a much better approximation. Nevertheless, any quantitatively 
reliable approximation should retain the one-gluon exchange in addition to the instanton 
interaction and the NJL-interaction (II . 11) . 



We emphasize that our approximations for the flow equation in the nonperturbative 
range remain relatively crude, even if in Sect. 3 we have done some checks in order to test 
the validity of some of the many approximations involved. Our results should therefore 
only be trusted on a semi-quantitative level. In our opinion, they are nevertheless a good 
guide for the qualitative characteristics of the effective four-quark interactions at scales 
between a few hundred MeV and 1 GeV. 
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FIGURE CAPTIONS 



Fig. 1. One-loop diagrams for the flow equation fl2.31| ) for the four-fermion (1PI) coupling 
A(j. Five more diagrams similar to c) with gluon lines attached at different quark 
lines are not shown. 



Fig. 2. Scale dependence of the quartic four-quark coupling. 

Continuous line: X a (t = 0) as a function of k with the initial condition as in Eqs. 
fl3.1| ), ( |3.7| ). Dotted line: X a (t = 0) as a function of k with the initial condition 
\ a (t)\ k=kc = [ansatz (A)]. Dot- dashed line: the same as for the continuous line, 
but using everywhere (both in the flow equation and in the initial condition) the 
momentum-independent coupling constant a [ansatz (B)]. Dashed line: the one- 
gluon-exchange contribution X£ 9 \t = 0), defined by Eq. ( |2.41| ), as a function of 
k. Long-dashed line: the instanton contribution X^ sU defined by Eq. (J4.3|), as a 
function of k. 



Fig. 3. Different approximations for X a in perturbation theory. 

Continuous line: the same as in Fig. 2. Dotted line: \&\t = 0) versus k, momentum- 
dependent a [ansatz (A)]. Dashedline: X^(t = 0) versus k, momentum-independent 
a [ansatz (B)]. Dot- dashed line: X^'(t = 0) versus k, constant a [ansatz (C)]. 

Fig. 4. Dependence of the four-quark interaction on the transferred momentum t. 

We show X a (t) (continuous line) and the one-gluon-exchange contribution A™(i) 
(dashed line) as functions of t at the fixed value k = 0.42 ■ Aqcd/ A^ GeV. 

Fig. 5. Quartic coupling with spontaneous color symmetry breaking. 

Continuous and dashed lines: the same as in Fig. 2. Lower (A) and upper (B) 
dotted lines: X a (t = 0) versus k using the running gluon mass term m\ given by Eq. 
(p|) (A) and Eq. Q ( B )- Lower ( A ) and u PP er ( B ) dot-dashed lines: A^(t = 0) 
versus k using the running gluon mass term m\ given by Eq. ( |5.2| ) (A) and Eq. 
(PD (B). 
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Fig. 6. Momentum dependence of quartic couplings with spontaneous color breaking. 

Continuous and dashed lines: the same as in Fig. 4. Lower (A) and upper (B) dotted 
lines: X a (t) versus t using the running gluon mass term m\ given by Eq. ( |5.2| ) (A) 
and Eq. ( |5.3| ) (B). Lower (A) and upper (B) dot-dashed lines: \£ 9 \t) versus t using 
the running gluon mass term ml given by Eq. ( |5.2| ) (A) and Eq. ( |5.3| ) (B). 
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